BCS-BEC crossover with unequal mass fermions 
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We investigate the crossover from BCS pairing to molecular BEC in an atomic gas with two 
fermion species with masses m-f 7^ m| tuned through a Feshbach resonance. We present results 
for the T = equation of state as a function of the scattering length including the effects of 
Gaussian fluctuations about the mean field ground state. We compute the ground state energy as a 
function of m^-/m^ at unitarity and find excellent agreement with the quantum Monte Carlo result 
for m-f /mi — 6.67 for a *°K-^Li mixture. We show that the dimer scattering length in the BEC 
limit as a function of m|/m| compares well with exact four-body results of Petrov et al. We also 
derive the condition for trapping frequencies to obtain an unpolarized gas in a harmonic trap. 



The crossover [H [1, Q of collective BCS pairing to the 
BEC of tightly bound molecules is of great interest in di- 
verse areas of physics ranging from condensed matter to 
nuclear and high-energy physics. The unitary scattering 
regime which lies in the middle of the crossover in three 
dimensions permits us to gain a deeper understanding of 
superfluidity in a strongly interacting Fermi system. Re- 
cent advances in atomic cooling and trapping and the use 
of a Feshbach resonance to tune the interactions through 
unitarity have led to the experimental realization Q of 
the BCS-BEC crossover in ultracold atomic gases. 

In this letter we turn to the question of pairing and 
superfluidity in a mixture of two species of fermions with 
different masses. There are several motivations for un- 
dertaking such an investigation. First, from an exper- 
imental point of view, interspecies Feshbach resonances 
have already been observed in mixtures of fermionic ^Li 
and ^''K This provides a whole new system in which 
strong interactions in the BCS-BEC crossover can be 
studied. Second, in high energy studies of pairing, such as 
in color superconductivity, the pairing naturally occurs 
between fermions with different masses Third, from a 
general many-body physics point of view it is interesting 
to ask if the unequal "spin" t and J, mass, which breaks 
"time-reversal" symmetry in condensed matter systems, 
affects the pairing of |k, t) and | — k, J,) states. 

Here we present results for the T = equation of 
state for an unpolarized system for arbitrary scattering 
length and mass ratio. We generalize to the case of un- 
equal masses a functional integral formulation Q that 
goes beyond the mean field theory results reported re- 
cently [8| and incorporates the effects of quantum fluctu- 
ations. Specifically, our approach includes the energy 
of the zero point motion of the collective mode (i.e., 
the Goldstone mode in the broken symmetry superfluid 
state) as well as the effects of virtual scattering of two- 
particle excitations. 

Our results (Figs. 1,2,3) are significantly different from 
mean field theory and also from those of the 1/J\f ap- 
proximation but are in very good agreement with 
the few known exact results. At unitarity, our results 
compare well (Fig. 2) with quantum Monte Carlo results 



that are available for two values of the mass ratio 
and Li-K mixtures 



HI 



equal 

In the BEC limit we 



masses [1 

find that our dimer scattering length as a function of the 
mass ratio (Fig. 3) is in very good agreement with the 
exact four-body results [13] ■ In the BEC regime we also 
find approximately 90% of the the Lee- Yang-Huang cor- 
rection caused by quantum depletion of the condensate. 
The Hamiltonian density for the system is given by 



2m^ 



where we use "spin" a 
species with masses 



(1) 

to label the two fermion 
> mi. The short-range attrac- 
tion |13| with ultraviolet cutoff A has a strength 5(A) 
that is related to the s-wave scattering length via [3| 
rn/Anas = —l/g{A) + X)|k|<A "^/^^i '^^^'^ defined be- 
low. We set h = 1 and work in a box of unit volume. 

We have two chemical potentials fi^ to obtain the den- 
sities rto- — n/2, even in the unpolarized case; see be- 
low. It is convenient to use /i — (fii + ^t)/2 and h — 
(/i| — AiT)/2 ™ place of the /i^, and also to define the 
reduced mass m/2 and < 7 < 1 given by 



m mimi 
2" ~ 



TO| + mi 



and 7 



m| — mi 
mi +mi' 



(2) 



We briefiy describe our formalism Q to intro- 
duce notation and highlight the changes arising 
from mi 7^ m^. The partition function Z(T,iJ,f^) 
is written as a functional integral over Grassmann 
fields. Introducing a Hubbard- Stratonovich field A(a;) 
which couples to ■ipi{x)ip i{x) and integrating out the 
fermions we obtain Z = J DADA*e~^ with the action 
S = /dx(|A(x)|V5-TrlnG-i[A(a;)]). The inte- 
gration is over x — (x, r) where the imaginary time 
< T < /3 = 1/T. The inverse Green's function is = 
[{~dr -h- 7VV2TO)i + (VV2m + ^i)f3 + Af + + A*f- 
xS{x — x'), and the trace acts on Nambu-Gorkov space. 

Mean field (MF) theory corresponds to a uniform, 
static saddle point Aq for Z. We find the gap equa- 
tion SSq/SAq = with 5*0 = S[Ao] = /3Ag/.g - 



E 



M.„TrlnGo-^(fc). 



Here '-{k) = {ik,, - h + je)l - 



f-Tj + AqTi with fermionic Matsubara frequencies ikn = 
i{2n + l)7r//3. To simplify notation we omit k labels on 
all energies and write them as e = |kp/2m and ^ = e — ^. 

The fermionic excitation energies, determined from the 
poles of Go (k, (jj + zO"*" ) , are 



Ei,2^ET{ie~h) with £;= + A2)1/2. 



(3) 



The thermodynamic potential = — In Z//3 in MF the- 
ory is given by Qq = So/P = Al/g + EkK " + 
EiQi^Ei) + E2e{-E2)] at T = 0. 

Precisely at T — 0, given /j,, Aq and 7, there is a range 
of h values for which we obtain an unpolarized solution 
n| = with E12 > 0. This range is max(7e ^ E) = 
< h < = min(7e-|-i^), where the max/min are 
for e > 0. It is easy to see that A^^^^ — ^fi : 
for fi/Ao > ±r while 

±r, 



^v/r^Ao 

- ±v//i' + A^ for fi/Ao < 
In the limit T ^ there 
+ A(~)]/2; deviations 



where F = 7/ a/1 — 7^ 
is a unique solution [l^ /i = 
away from this value lead to an exponentially small po- 
larization at low T within the interval. 

The expression for Sq simplifies greatly with £^1^2 > 
and we obtain the gap equation 



(4) 



and MF number equation n 



^^(1 — ^/i?). Thus the properties of the m-^ 7^ 
system at the MF level are completely equivalent to 
those of an equal mass system with chemical potential 
fi = [fii + fi^)/2 and m twice the reduced mass ([2]). This 
was already noted in ref . 0] ; here we go beyond their MF 
analysis by including the effect of Gaussian fluctuations 
on the equation of state. 

We write A(x) = Aq -I- where the -q are (x, r) de- 
pendent fluctuations about the saddle point Aq. Fourier 
transforming to momentum q and bosonic frequencies 
iqi = i27r€//? and expanding to quadratic order in 
rj we get 5*^50 + Sg, where Sg = ^ I^q,^?, '/^M?7 
with ~ [rj* {q),ri{~q)) The inverse fluctuation prop- 
agator M is given by Mii(g) = M.22{~q) — I/5 + 
r'Ek^fe G°22ik)GUk + q) and M,2{q) = M2i{q) = 
Ek .fe„ G?2(fe)G?2(fc + q). At T = we get 



Mn(g) 
Mi2(g) 



V V 



iqi-Ei-E'2 iqe + Ei + E'^ 



iqe + El 



E!2 



El — E2 



(5) 



Here we use standard BCS notation v"^ = 1 — = 
i (1 — f/i?) together with the abbreviations u — Uk,u' = 
Uk+q, etc. The only difference in from the equal 
mass case is the presence of i?i + i?2 = + -E' — 7 (e' — e) . 
Thus, even at the Gaussian level, h does not enter the 
calculation so long as the system is unpolarized. 
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FIG. 1; Chemical potential /i normalized by ep = (eft + 
epi)/2 as a function of the interaction — l/(fcFas) with 
nil /mi = 6.67 for ^"K and ^Li. The solid (red) line includes 
Gaussian fluctuations, while the dashed (black) line is the 
mean field prediction (independent of mi /mi). 



The partition function is now approximated by Z ~ 
exp(— S'o) / Dr]Dr]^ exjp{—Sg). As explained in ref. 0, 
Il5j the saddle point equation still retains its earlier form 
SSq/SAq = 0, but the thermodynamic potential 17 at T = 
now includes the contributions of the energy of zero 
point motion of the collective excitations (the Goldstone 
mode of the broken U{1) symmetry in the superfluid) and 
the virtual scattering of gapped fermionic quasiparticles. 
A careful consideration of convergence factors [t} leads 



to SI ~ r^o -f iln 



with 



= ^ E ln[Mn(g)DetM(g)/M22(q)]. 



(6) 



The Gaussian contribution leads to a lowering of the 
ground state energy and a modified equation of state 
which we calculate as n = —dfl/dfj, at T = 0. Given 
the dependence of Aq on /i, we thus find 



n = —dVlo/dfi — 9rig[/i, Ao(^)]/9^. 



(7) 



In Fig. 1 we show the average chemical potential /i as a 
function of the scattering length for the experimentally 
relevant case of a *°K-^Li mixture. 

To understand in greater detail how unequal masses 
alter the results from the well studied equal mass case, 
we next concentrate on two regions: unitarity and the 
BEC limit. It is also useful to compare our results with 
those of a closely related approach, where one general- 
izes the Hamiltonian ^ to M flavors of fermions of each 
species, interacting with an Sp(2A/')-symmetric poten- 
tial [9] . Then MF theory becomes exact in the Af ^ 00 
limit and the Gaussian correction is of order l/A/" so 
that n/J\f ^ flo + VLg/N + 0[1/N^). Thus, we write 
A* = /^o + Ml /A/" where /io is the MF value obtained from 
n = —dilo/du, and fii = ~{dVLg/dfi)/ {d'^VLo/ dfx^) eval- 
uated at /iQ [l6[. In order to obtain results relevant to 
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FIG. 2: The chemical potential at unitarity as a function 
of the mass ratio. At unitarity, the ground state energy E 
and /i are related by universality: t-i/ep = E/ {SNIf /5) — 
(1 + /3) where If = {ep^ + epi)/2. Our calculation including 
Gaussian fluctuations is the solid red line. The vertical line 
marks /mi = 6.67 for a ^"K/^Li mixture. For comparison, 
the two triangles are Quantum Monte Carlo results the 
upper dashed line is the mean field result and the lower dotted 
line that of the 1/N expansion. 



the original Hamiltonian M — 1 \s set at the end of the 
calculation. In contrast, in our Gaussian theory we work 
directly with the Hamiltonian ([T]), keeping A/" = 1 from 
the outset, and calculate corrections to various quantities 
(like the chemical potential) without assuming that the 
changes are small. 

At unitarity 0^=00 and the interactions have no scale; 
the chemical potential is the only energy available 17 1. 
The saddle point condition implies ^ = 0.860 Aq inde- 



pendent of mass ratio and the thermodynamic potential 
has the universal scaling form 



■J^(m|/m|)/i^/2 (2^)3/2 



(8) 



1 for 



where the function J- is normalized so that 
non-interacting fermions with equal masses. 

In the MF approximation = 2.2032, independent of 
ra^/mi, which yields the equation of state //q = ^f/^o^'^ 
where e_F = kp/{2m) = (e^-f +eFj,)/2. The Gaussian con- 
tribution J-g has a non-trivial, though weak, dependence 
on the mass ratio obtained by evaluating ^ numerically. 
This yields n = ep/l^o + •^g("^T/"^i)]^^'^- This result is 
plotted in solid line in Fig. 2, with the dashed line show- 
ing the mean- field result. The dotted line corresponds 
to the l/A/" result ^ — [I — 2Tg{m-\/mi)/iN!FQ\^Q with 
J\f = 1. In Fig. 2 we also show the result of a recent 
Quantum Monte Carlo calculation [ll[ (triangles) which 
agrees fairly well with our approximation. 

Let us next discuss the BEG limit (a^ ^0+) where we 
obtain a condensate of tightly bound diatomic molecules 



of binding energy Eb = 
and density ub = n/2. 



= l/{mai), mass vtlb = rn^ + rni 
These bosonic dimers are weakly 



6- 
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FIG. 3: Molecular scattering length in the BEC limit as a 
function of mass ratio for different calculations. Our calcula- 
tion including Gaussian fluctuations is the solid red line; the 
exact 4-body calculation [l3| is shown as the dot-dash line. 
The mean field value, corresponding to the Born approxima- 
tion for dimer scattering, is the upper dashed line. 



interacting with a (dimer) scattering length a^, which 
we now compute. We find it convenient to define di- 
mensionless variables x and y as \^\ = (1 — y)/{2mal), 
1 and (2ma^Ao)^ = a; ^ 1 in the 
We can thus expand to quadratic order 
[—xy + x^ /16]/{32nmal). The gap equa- 



with < y < 1 and (2ma^Ao) 
BEC hmit 
SQ{x,y)/(3 
tion dSo/dx = then yields x — 8y. 

At the MF level f7o(/u) = Soix = 8y,y)/l3 so that the 
number equation n = — (2ma'^)dflo / dy yields y = 2Trna^. 
We find /i ~ —l/{2ma'^) + Trnas/m and on identifying the 
second term with /zb/2 = 2TrnBaB /itt-b we obtain the 
MF result [1] a^^/us — {m^ -|-TO|)/to. This corresponds 
to the Born approximation for dimer scattering and is 
shown as a dashed line in Fig. 3. 

The Gaussian contribution © in the BEC limit can 
be written as fig = — x'^g{m^ /mi)/ [512Trma^] where the 
function g{mi/mi) can be numerically calculated by the 
method of ref. 0] (Appendix E). Using x = 8y,we write 
ri(^) = Qo+rig = —{l+g)y'^ / (Smna^g). Solving the num- 
ber equation n = — (2ma^)dfl / dy for ^ we determine the 
dimer scattering length (as explained above) to find [l^ 
as/as = ("^T + + 5("^T/"^i))"^]- This is plotted 

as the solid line in Fig. 3 and compares very well with the 
exact result from a 4-body calculation of the scattering 
between dimers 



12| 



We use our numerical results to find the next correction 
in the BEC limit, by fitting fl^n) at small y to to the 
form {Ay'^ + By'^/'^ + . . .)/{ma^^). We thus find the T = 
equation of state of the Lee- Yang-Huang form 

t^B = + 32C(nsas)'/V(3V^) + ■ • ■] (9) 

TUB 

arising from the quantum depletion of the condensate. 
From our numerics we find C — 0.90 ±0.05 over the mass 



4 



ratio range 1 < m|/mx < 13.6. The exact result for 
the dilute Bose gas is C = 1. 

Finally we analyze the effect of a harmonic trap within 
the Thomas- Fermi local density approximation (LDA); 
for simplicity we discuss spherical traps. The two species 
will in general see different trapping potentials with fre- 
quencies oja-, SO that within LDA /icr(r) = Jl^—ma-uj'^r'^ /2. 
Let us focus on unitarity where Ao(r)//i(r) = 1.162 = 6 
for all r. To ensure zero polarization at T = we must 
satisfy h{r)/fi{r) — K. Using the results derived below 
eq. (3), X = 7 for 7/^1 — 7-^ < 1/^ corresponding to 
m^/mi < 4.74, while X = [7 - J^/l - 7^ + Vl + S'^]/2 
for TO| / mi > 4.74. To maintain this ratio of we must 
choose 



(10) 



Thus, LOi/uji = mi /mi for m^/mi < 4.74 while for a 
Li-K mixture, where K = 0.745, we get coi/coi — 0.148. 
Deviations from this frequency ratio would lead to re- 
gions of spin imbalance in the trap. 

We conclude with a brief discussion of some general 
issues. First, note that mass imbalance does not im- 
pact pairing of "time-reversed" states. For the unpo- 
larized case the Fermi spheres match in k-space even if 
the ^'s are different. Unlike in many condensed matter 
systems where time-reversal breaking perturbations [20*1 
flip the spin of the fermions, atomic species retain their 
identity ("spin") in collisions. However, effects beyond 
pairing (which is what we focused on) may well desta- 
bilize the system. In the BEC limit, it is known that 
for mi /mi > 13.6 the presence of Efimov states makes 
the Feshbach molecules unstable to collapse by three- 
body recombination [l^ . Finally, we may ask why our 
approach seems to work so well when compared with 
available exact results (at least in the stable mass ra- 
tio regime). The answer seems to lie in the fact that 
in a homogeneous dilute gas (with kp^ ^ range of po- 
tential) particle-hole channel contributions are negligible 
compared to the particle-particle channel ones that our 
approach is designed to capture. 
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